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The ding-a-ling model is a kind of half lattice and half hard-point-gas (HPG) model. The original 
ding-a-ling model proposed by Casati et. al does not conserve total momentum and has been found 
to exhibit normal heat conduction behavior. Recently, a modified ding-a-ling model which conserves 
total momentum has been studied and normal heat conduction has also been claimed. In this work, 
we propose a full lattice ding-a-ling model without hard point collisions where total momentum is 
also conserved. We investigate the heat conduction and energy diffusion of this full lattice ding-a- 
ling model with three different nonlinear inter-particle potential forms. For symmetrical potential 
lattices, the thermal conductivities diverges with lattice length and their energy diffusions are su- 
perdiffusive signaturing anomalous heat conduction. For asymmetrical potential lattices, although 
the thermal conductivity seems to converge as the length increases, the energy diffusion is definitely 
deviating from normal diffusion behavior indicating anomalous heat conduction as well. No normal 
heat conduction behavior can be found for the full lattice ding-a-ling model. 

PACS numbers: 05.60.-k,44.10.+i,05.45.-a 


I. INTRODUCTION 


It has long been believed that there should be profound 
connection between heat conduction and energy diffusion 
since both of them describe the same dynamic process of 
energy transport 'lj. For normal and ballistic heat con¬ 
duction, it is well established that their energy diffusions 
should also be normal and ballistic, respectively. What 
is nontrivial and really interesting is the connection be¬ 
tween anomalous heat conduction and anomalous energy 
diffusion for low dimensional systems. For anomalous 
heat conduction, the thermal conductivity diverges with 
the system length as k oc N a with 0 < a < 1 [ i - l 3 |. For 
anomalous energy diffusion, the spreading of the Mean 
Square Displacement (MSD) of energy follows the be¬ 
havior as (A x 2 (t)) E oc with 1 < /3 < 2 . In 

the early pioneering studies, it has been found that the 
connection formula between heat conduction and energy 
diffusion is a = (3— 1 . Only until recently, a rigor¬ 

ous connection theory between heat conduction and en¬ 
ergy diffusion has been established within the framework 
of linear response theory [ 40 .]. As a byproduct, the above 
mentioned connection formula comes out quite naturally. 


According to the connection theory [4(|, the study of 
the energy diffusion is equivalent to the study of heat con¬ 
duction in the sense of calculating thermal conductivity. 
As has been pointed out by Ref. [U| , the temporo-spatial 
distribution of energy fluctuation correlation function 
gives more information than the heat current-current cor¬ 
relation function used in Green-Kubo formula. This jus¬ 
tifies the advantage of using energy diffusion method to 
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study the problem of heat conduction. In previous works, 
it has been found that the energy fluctuation correla¬ 
tion function is Gaussian for systems with normal heat 
conduction [Hj]. More interestingly, the energy fluctua¬ 
tion correlation function for systems with anomalous heat 
conduction is a Levy-walk like distribution [30l432l . [H- 
[37l Ell [42 ]. For a comprehensive understanding of the 
Levy wal ks, p lease refer to the following excellent refer¬ 
ences [H, |44| . In general, the ID lattice system exhibits 
normal heat conduction behavior if the total momentum 
is not conserved and anomalous heat conduction behav¬ 
ior if the total momentum is conserved EHH There is 
however one exception for the ID coupled rotator model 
which conserves total momentum but shows normal heat 
conduction behavior [H, 46] . Most recently, it has been 
confirmed that the ID coupled rotator model possesses 
normal momentum diffusion as well as normal energy dif¬ 
fusion m m, m. 

The first ID ding-a-ling model proposed to study the 
heat conduction problem was introduced by Casati et.al 
in 1984 |49]. This system consists of on-site harmonic os¬ 
cillators and free moving particles positioned alternately. 
It is a kind of half lattice and half HPG model where total 
momentum is not conserved. Therefore it is not a sur¬ 
prise that normal heat conduction has been found for this 
model [49[. Recently, a modified ding-a-ling model was 
proposed where the on-site harmonic oscillators are re¬ 
placed by inter-connected harmonic oscillators [50]. With 
this setup, the total momentum is conserved. It is then 
quite surprising that normal heat conduction has also 
been detected for this modified half lattice and half HPG 
ding-a-ling model [5Cj . 

In this paper, we will propose a full lattice ding-a- 
ling model which conserves the total momentum. Three 
different inter-particle potentials will be applied to this 
model. The heat conduction and energy diffusion be- 
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FIG. 1: The schematic setup for the full lattice ding-a-ling 
model conserving total momentum. The top particles are 
connected by harmonic springs. The bottom particles are 
connected with its two neighbors of top particles with sym¬ 
metrical and asymmetrical nonlinear inter-particle potentials. 


haviors will be studied for each potential form with two 
numerical methods. Firstly, the length dependent ther¬ 
mal conductivities k will be calculated using the non¬ 
equilibrium molecular dynamics simulations. Secondly, 
the energy fluctuation correlation function will be calcu¬ 
lated at different correlation times and the MSD of energy 
spreading will be presented to reveal the heat conduc¬ 
tion behavior for this momentum conserving full lattice 
ding-a-ling model. This paper will be organized as the 
followings: in Sec. II we will give the introduction of the 
full lattice ding-a-ling model with three different inter¬ 
particle potentials. The heat conduction and energy dif¬ 
fusion behaviors will be presented in Sec. III. In Sec. IV, 
we will give our summary for the study of heat conduc¬ 
tion as well as energy diffusion for our new proposed full 
lattice ding-a-ling model. 


II. ID MOMENTUM CONSERVING FULL 
LATTICE DING-A-LING MODELS 


The full lattice ding-a-ling model is depicted in Fig. 
[B The original half free moving particles shown as the 
bottom particles in Fig. [T|are now connected with its two 
top neighbors via the nonlinear inter-particle potentials. 
The Hamiltonian for this full lattice ding-a-ling model is 
the following: 


h = Y. h < = T. 


Pn.Pl, 1 , ,2 

— + + 2 Wt+1 “ 9iJ 

+v(<k b ,qi t ) + git + l)] 


( 1 ) 


where the index it and ib denote the particle index in 
the top and bottom respectively. The pt t and pt b are the 
momentum for top and bottom particles. The qt t and 
qi b are the displacements from equilibrium for top and 
bottom particles. 

The inter-particle potential V (qt t , ) describes the in¬ 
teraction between top and bottom particles which could 
be any nonlinear potential form. Here we consider three 



FIG. 2: (Color online) The schematic picture of the asymmet¬ 
rical half FPU-/3 form potential V ( qt b —qi t ) and V ( qt b —qt t + i). 


typical nonlinear potentials for V: (i) symmetric FPU-/3 
potential with 


V(q ibl q it ) 


(Qh-Qu) 2 fe b -9u) 4 
2 4’ 


(ii) symmetric coupled rotator potential with 


( 2 ) 


V{qi„,qi t ) = 1 - cos(q ib - q u ); (3) 

(iii) asymmetric half FPU-/3 potential depicted in Fig. [2] 
with 


V(q ib ,q it ) 

( (git ~ Qn) 2 , 

n 2 

V{q ib ,q it + i) 



-gn+1) 2 

2 


4 


Qi„ 


0 , qi b 


Qit < 0 

qu > 0 


( 4 ) 



o, qi b - gi t +i < o 


- gi t +1) 4 

4 


qi b - git+i > o 


The first two models are symmetrical since V ( qt b , qt t ) 
and V(qi b , ft t +i) posses the same form while the third 
model is asymmetric as V(qi b ,qi t ) and V(qi b . qi t +i) have 
different potential forms. Compared with the previously 
studied half lattice and half HPG ding-a-ling model, 
this full lattice ding-a-ling model enables us to simu¬ 
late the system dynamics with more accurate and consis¬ 
tent numerical algorithm such as fourth order symplectic 
method [5l| in studying the energy diffusion behaviors. 
In the following calculation, the dimensionless units have 
been applied. 


III. HEAT CONDUCTION AND ENERGY 
DIFFUSION 

In the study of heat conduction, we first apply the di¬ 
rect non-equilibrium molecular dynamics method. For 
simplicity and convenience, the fixed boundary condi¬ 
tions with qt t -o = <7j t= jv +1 = 0 will be used. The first 
and last atoms on the top chain will be contacted with 
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FIG. 3: (color online). Thermal conductivity k as the function 
of lattice size N for symmetric FPU-/3 potential. The circles 
are the numerical data and the straight line of k oc IV 0 ' 45 is 
guided for your eyes. The average temperature is chosen as 
T = 12.15 corresponding to energy density E = 10 used in 
equilibrium method below. Temperature bias is set as A = 
0.1 and N is the number of unit cells. 


two Langevin heat baths with temperature T r /r. In par¬ 
ticular, the equations of motions for these two atoms are 


Qu = F {Qn) - Mu + (*) (5) 


where F(qi t ) is the force which can be derived from 
Eq. Q and corresponding boundary conditions, and the 
Gaussian white noise £* t =i/jv(t) satisfy 

(£i t =i/jv(*)) = 0 

<4=i/jv(i)6 t =i/jv(0)) = 2A k B T L/R 5{i) (6) 

where A is the damping parameter which will be always 
set as unity and (■) denotes the ensemble average which is 
equivalent to the time average for the chaotic systems we 
considered here. The second order Verlet velocity method 
will be used to integrate the equations of motions. 

The unit cell contains one top particle it and one bot¬ 
tom particle if, indexed by i = it = ib- From continuity 
equation of energy for this unit cell, the heat flux can be 
derived as j* = -q it {q it - q it - 1 ) - qu d V{q ib ^i,q it ) / dq it . 
The thermal conductivity can be calculated as 


where J = {j t ) is the average heat flux in the stationary 
state independent on the index i and VT is the tempera¬ 
ture gradient. The temperature of the heat baths are set 
as T l / r = T{ 1 ± A) where T is the average temperature 
and A is the bias. 

In the study of energy diffusion, the equilibrium nu¬ 
merical method will be applied with the periodic bound¬ 
ary conditions of q H = <jN+i t and q ib = qw+u ■ The key 




FIG. 4: (color online), (a) Energy fluctuation correlation 
function Gfi(i,t) for model with symmetric FPU-/3 potential. 
The red, olive and blue solid lines represent the functions at 
t = 100, 200 and 300, respectively, (b) The MSD of the excess 
energy (Ax 2 (t)) E as the function of correlation time t. The 
two solid reference lines of t 1 ' 45 and t are guided for the eyes. 
It can be seen that the numerical data almost follows the 
t 1 ' 45 behavior which is the similar anomalous energy diffusion 
behavior previously found for ID FPU-/3 lattice. The energy 
density E = 10 and the lattice length N = 1001. 


information needed to be calculated is the energy fluctu¬ 
ation correlation function CE{i,t) defined as [35| 


(AHj(t)AH o (0)) 

(Affo(O)Affo(O)} 


( 8 ) 


where A Hi{t) = Hi{t) — {Hi). The local energy Hi in¬ 
cludes the energy of top particle i t and bottom particle 
ib with i = i t = ib- The unit cell index i runs from 
— {N — l)/2 to {N — l)/2 so there are N top and bottom 
particles each. With this indexing for odd N, the central 
particles have the index it = ib — 0 which turns out to 
be convenient for numerical simulations. The system will 
have 2 N particles in total. Note that CE{i,t = 0) = <5^0 
is the initial distribution for energy fluctuation correla¬ 
tion function. If one consider a small initial excess energy 
perturbation with a special delta form, the excess energy 
distribution pE{i,t) equals to the energy fluctuation cor¬ 
relation function as pE{i,t) = CE(i,t), which describes 
the actual time evolution of the initial excess energy along 
the lattice [iol ]. 

The MSD of the excess energy distribution (Arr 2 (t)} £ 
can be defined as [H, S3 : 


(A x 2 {t)) E = Y^i 2 P E {i,t) = ^2i 2 C E {i,t) (9) 

i i 

According to the connection theory, the second derivative 
of the MSD of the excess energy distribution (Ai 2 (t)) £; 
is connected with the autocorrelation function of total 
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FIG. 5: (color online). Thermal conductivity k as the function 
of lattice size N for symmetric rotator potential. The circles 
are the numerical data and the straight line of k oc IV 0 ' 48 is 
guided for your eyes. The average temperature is chosen as 
T = 0.91 corresponding to the energy density E = 1 used in 
equilibrium method below. Temperature bias is set as A = 
0.2 and N is the number of unit cells. 


heat flux Cjj(t) (40l | : 

d 2 (Ax 2 (t)) E _ 2 Cjj(t) 

dt 2 k B T 2 c 1 ' 

where c is the volumetric specific heat capacity and k B is 
the Boltzmann constant. The autocorrelation function of 
total heat flux Cjj (t) is the central quantity which enters 
the Green-Kubo formula for the calculation of thermal 
conductivity Si: 

-■i&L cjAt)dt (u) 

The connection theory between heat conduction and 
energy diffusion tells us the following [40]]: if the energy 
fluctuation correlation function is Gaussian in 

the asymptotic time limit, the MSD of the energy distri¬ 
bution (Ax 2 (t)} E will be linearly proportional to time 
t as (Ax 2 (t)) E oc t. This linear time dependence of 
(Ax 2 (t)) E eventually gives rise to a finite thermal con¬ 
ductivity k indicating a normal heat conduction behav¬ 
ior. This has been already observed for </> 4 lattice [35} and 
coupled rotator lattice [4T( where both of the systems ex¬ 
hibit normal heat conduction [52], [53} ■ If the energy fluc¬ 
tuation correlation function C B (i,t) is Levy-walk like in 
the asymptotic time limit, the MSD will grow faster than 
linear time dependency as (A x 2 [t)) E oc t 13 with /? > 1. 
This is the case for ID momentum conserving FPU-/? 
lattice, amended rotator lattice and Lennard-Jones lat¬ 
tice [llllil]]. From the connection theory [40f. the thermal 
conductivity will diverge as k oc N a with a = (3 — 1 ex¬ 
hibiting anomalous heat conduction behavior. 

In order to study the heat transport behavior for the 
full lattice ding-a-ling models with three different non- 



FIG. 6: (color online), (a) Energy fluctuation correlation 
function Gfi(i,t) for model with symmetric rotator potential. 
The red, olive and blue solid lines represent the functions at 
t = 300, 400 and 500, respectively, (b) The MSD of the 
excess energy (A x 2 (t)) E as the function of correlation time t. 
The two solid reference lines of t 1 ' 48 and t are guided for the 
eyes. It can be seen that the numerical data also follows the 
t 1 ' 48 behavior which is the similar anomalous energy diffusion 
behavior previously found for ID FPU-/3 lattice. The energy 
density E = 1 and the lattice length N = 1001. 


linear inter-particle potentials, we will use the two nu¬ 
merical methods to investigate the heat conduction as 
well as the energy diffusion behaviors. For the non- 
equilibrium molecular dynamics simulations, the length 
dependent thermal conductivities n will be calculated di¬ 
rectly. For the equilibrium simulations, we will calcu¬ 
late the temporo-spatial distributions of energy fluctua¬ 
tion correlation functions C B (i,t ) at different correlation 
times. The resulted MSD of the excess energy (Ax 2 (f)) £ 
as the function of correlation time t will be used to an¬ 
alyze the actual behavior of energy diffusion for each 
model. 


A. Symmetric FPU-/3 potential 

We first consider the full lattice ding-a-ling model with 
symmetric FPU-/? potential. The length dependent ther¬ 
mal conductivities k are plotted in Fig. [3] with average 
temperature T = 12. The k shows a good power-law 
length dependence as k oc N 0A5 . The heat conduction 
should be anomalous similar to the FPU-/? lattice. Simu¬ 
lations with other average temperatures show almost the 
same results. 

In Fig. Eta), the energy fluctuation correlation func¬ 
tions Ce(i, t) with symmetric FPU-/? potential have been 
plotted at different correlation times t = 100, 200 and 
300. Similar to the most studied ID FPU-/3 lattice [35} . 
the correlation function C B {i, t) here is also like the Levy- 
walk distribution which is characterized by one central 
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FIG. 7: (color online). Thermal conductivity k as the function 
of lattice size N for asymmetric FPU-/3 potential. The circles 
are the numerical data and the thermal conductivity k seems 
to converge in the larger N region. The average temperature 
is chosen as T = 1.0 corresponding to the energy density 
E = 1 used in equilibrium method below. Temperature bias 
is set as A = 0.2 and N is the number of unit cells. 

peak and two side peaks. The side peaks moves towards 
outside with a constant sound velocity. To determine the 
diffusion behavior, the time dependence of MSD of the 
excess energy (A x 2 (t)) E needed to be calculated. The 
MSD of the excess energy (A x 2 (t)) E as the function of 
correlation time t has been plotted in Fig. 01b). As can 
be seen from the figure, the MSD of the excess energy 
exhibits superdiffusion behavior as (A x 2 (t)) E oc t 1A5 . 
Our results are consistent with the previously found dif¬ 
fusion exponent /3 = 1.40 for the ID FPU-/3 lattice [35| . 
According to the connection theory between energy dif¬ 
fusion and heat conduction|40|. superdiffusion of energy 
of /3 = 1.45 will imply an anomalous heat conduction 
with k oc 7V'“ =0 - 45 for the considered full lattice ding- 
a-ling model with symmetric FPU-/3 potential. This 
momentum-conserving ding-a-ling model exhibits almost 
the same energy diffusion behavior as the ID momentum- 
conserving FPU-/3 lattice. We have also done the simula¬ 
tions for this model at other energy densities and similar 
results are found. 


B. Symmetric rotator potential 

We then consider the full lattice ding-a-ling model with 
symmetric rotator potential. The length dependent ther¬ 
mal conductivities k are plotted in Fig. [5] with average 
temperature T = 0.9. The thermal conductivity k still 
shows a good power-law length dependence as n oc N 0A8 . 
The heat conduction is anomalous as previous case. Sim¬ 
ulations with other average temperatures show almost 
the same results. 

In Fig. Eta), we plot the energy fluctuation corre¬ 
lation functions C E {i,t) at different correlation times 
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FIG. 8: (color online), (a) Energy fluctuation correlation 
function CE(i,t) for model with asymmetric FPU-/? poten¬ 
tial. The red, olive and blue solid lines represent the func¬ 
tions at t = 200, 300 and 400, respectively, (b) The MSD 
of the excess energy (Ax 2 ^)} E as the function of correlation 
time t. The two solid reference lines of t 1 ' 80 and t are guided 
for the eyes. It can be seen that the numerical data follows 
the l 1 ' 80 behavior which is faster than the anomalous energy 
diffusion behavior previously found for ID FPU-/3 lattice, (c) 
The decay of the central peak C E (i = 0, t). It is obvious that 
the decay of the central peak as C E (i = 0, t) oc f -1 is faster 
than the decay for normal diffusion case which should be pro¬ 
portional to f -1 / 2 . The energy density E — 1 and the lattice 
length N = 1001. 


t = 300, 400 and 500 for the ding-a-ling model with 
symmetric rotator potential. The clear characteristics 
of Levy-walk distribution with one central peak and two 
side peaks appear. Compared with the ding-a-ling lattice 
with FPU-/3 potential, the side peaks here are relatively 
small. But the time dependence of the MSD of the excess 
energy still follows the same time dependent behavior as 
(A x 2 {t)) E oc t 1AS as can be seen from Fig. Etb). Ac¬ 
cording to the connection theory, the heat conduction for 
the ding-a-ling lattice with rotator potential is anomalous 
with n oc A 0 48 which is the case in Fig. El 

It is interesting that although the ID rotator lattice 
exhibits normal ener gy diffusion and normal heat con¬ 
duction behavior [4j], 07], HH , the ding-a-ling model with 
rotator potential between top and bottom atoms displays 
the similar super energy diffusion and anomalous heat 
conduction behavior as the ID FPU-/3 lattice [35j . 
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C. Asymmetric FPU-/T potential 

To further mimic the half lattice and half HPG ding- 
a-ling model where the free particle can only collide with 
the particle in its moving direction, we introduce the 
ding-a-ling lattice with asymmetric FPU-/3 potential as 
can be seen in Fig. [2] The thermal conductivities k are 
plotted as the function of length TV in Fig. [7] The k first 
increases with length TV and then seems to saturate to a 
constant value after TV > 1000. It looks like that a = 0 
and the heat conduction is normal in this asymmetric 
case. 

However, if we look at the energy fluctuation corre¬ 
lation functions CE(i,t) at different correlation times 
t = 200, 300 and 400 plotted in Fig. (HJa), the signa¬ 
tures of Levy-walk distribution with one central peak and 
two side peaks have also been found. The only differ¬ 
ences are the negative tips next to the side peaks which 
has already been found for other ID asymmetric lattices. 
The time dependence of the MSD of the excess energy 
has been plotted in Fig. [Bib) and a superdiffusion with 
(Aa; 2 ^))^ oc t 1 11 ’ 80 has been obtained. The relation be¬ 
tween a and (3 here deviates from the connection the¬ 
ory as noticed in recent work [54]]. It might be that the 
asymptotic length is very large for asymmetric lattices. 

We then further plot the decay of the central peak 
Ce{i = 0 ,t) in Fig. (He). It is observed that the decay 
rate is proportional to f -1 which is also much faster than 
the decay of the central peak for a normal diffusion which 
should be proportional to f -1 / 2 . Therefore, it will be 
insufficient to judge the heat conduction behavior only 
via the non-equilibrium molecular dynamics simulation 
for asymmetric potential lattices. One should rely on the 
energy diffusion behavior as well. 


IV. SUMMARY 

In summary, we have systematically studied the 
heat conduction and energy diffusion in momentum- 
conserving ID full lattice ding-a-ling model with sym¬ 
metric and asymmetric potentials. For symmetric lat¬ 
tices, the thermal conductivities diverge as the lattice 
length with a power-law dependence displaying obvious 
anomalous behaviors. The super energy diffusions of 
(A x 2 (t)) E oc with /? > 1 have also been found. For 
asymmetric lattices, the thermal conductivity seems to 
converge as the lattice length increases. However, the 
energy diffusion behavior is definitely superdiffusive in¬ 
dicating anomalous heat conduction behavior also. As a 
result, the heat conduction for the full lattice ding-a-ling 
models should be anomalous as that in ID FPU-/3 lat¬ 
tice @. This is in contrast to the momentum-conserving 
half lattice and half HPG ding-a-ling model where normal 
heat conduction has been obtained [50]. The underlying 
physical mechanism behind this discrepancy is an open 
issue and deserves further investigation in the future. 
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